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We study the ground state properties of bosons in a tilted double- well system. We use fidelity susceptibility to 
identify the possible ground state transitions under different tilt values. For a very small tilt (for example 10 -10 ), 
two transitions are found. For a moderate tilt (for example 10~ 3 ), only one transition is found. For a large tilt 
(for example 10 _1 ), no transition is found. We explain this by analyzing the spectrum of the ground state. The 
quantum discord and total correlation of the ground state under different tilts are also calculated to indicate those 
transitions. In the transition region, both quantities have peaks decaying exponentially with particle number TV. 
This means for a finite-size system the transition region cannot be explained by the mean-field theory, but in the 
large-TV limit it can be. 

PACS numbers: 67.85.-d, 03.75.Lm, 03.67 .Mn, 03.75.Gg 



I. INTRODUCTION 

The many -body quantum states in ultracold gases have been 
studied with a high interest, because there are many parame- 
ters that can be adjusted in experiments to control the static 
state as well as the dynamics of the system. As a paradigm 
model, Bose-Einstein condensate (BEC) in a double-well sys- 
tem provides a useful setup to tackle the properties of quantum 
systems. By loading ultracold atoms in double wells, one can 
study fundamental quantum mechanical effects and many im- 
portant quantum many-body phenomena, for example inter- 
ferometry [1], quantum information processing [2], quantum 
phase transition [3], quantum superposition state [4], Joseph- 
son oscillations and nonlinear self-trapping form of dynamics 
[5]. 

Recently, it was found that in some ultracold gas systems 
(such as rotating BEC and BEC in double wells with very 
small tilt), the ground state transition cannot be described by 
the mean-field theory, although the ground states before and 
after transition are consistent with the mean-field description 
very well [6-10]. This is because in the transition region, the 
ground state is no longer a product of single-particle states 
but a strongly-correlated entangled state. How to character- 
ize such states is under intensive study. Moreover, in Refs. 
[6-10], the system sizes under study are not too large. The 
property of quantum correlation in the transition region for 
larger system is an interesting problem. 

The appearance of quantum correlation in the ground state 
in the transition region makes it reasonable to use some tools 
borrowed from quantum information theory to investigate the 
transition. In this article, our goal is to use fidelity suscepti- 
bility and quantum correlation to study the ground state tran- 
sition of BEC in double wells with an arbitrary tilt. In Sec. II, 
we introduce the model and give the prediction of the ground 
state transition according to the semiclassical picture. In Sec. 
Ill, we calculate the fidelity susceptibility which can precisely 



'Electronic address: hlguophys@iphy.ac.cn 



locate the critical point of a possibly unknown quantum tran- 
sition [11]. Different behaviors of the fidelity susceptibility 
under different tilts are found. In Sec. IV, enlightened by the 
fact that entanglement can also show a rather interesting be- 
havior at the critical point of a quantum transition [12], we 
calculate two quantities of quantum correlation: total corre- 
lation and quantum discord, both of which can exhibit signa- 
tures of the quantum transitions [13]. Moreover, quantum dis- 
cord can appear even when entanglement is absent [14-18] so 
it is a more suitable quantity than entanglement to character- 
ize the quantumness of the correlation. We find that both total 
correlation and quantum discord are nonzero in the transition 
region, but their values decrease with the particle number. By 
doing a finite-size analysis, an exponential decay of them with 
the particle number is found. This means although for a small 
system the transition is dominated by the quantum correlation, 
for a very large system no quantum correlation exists during 
the transition. A brief summary is given in Sec. V. 



H. MODEL 

The single level Bose-Hubbard Hamiltonian for TV atoms 
in a double-well system can be written as 

H = - J(a^a R + a R a L )-U[h L (h L -l) + h R (h R -l)] 
- V (h L - h R ), (1) 

where a\ (di) creates (annihilates) a boson in the i-th well 
(i = L, R), hi = a\cii, J is the tunneling energy and U is the 
on-site interaction [a positive (negative) U corresponds to at- 
tractive (repulsive) atom-atom interaction]. Vq is the tilt which 
can destroy the left-right symmetry and is non-zero in real ex- 
periments. We set J — 1 for convenience and only consider 
U > in this paper. The above Hamiltonian can be diag- 
onalized in the (TV + 1) -dimensional Fock space spanned by 
\ n L, n R = N—ni). The dynamics of the system is controlled 
by the parameter A = NU / J. As A passes from the weak re- 
gion to the fermionization limit, the dynamics of these atoms 
which are initially prepared mostly in one well, will change 
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from Josephson oscillation (simply tunneling back and forth 
between two potential wells) to self-trapping above a critical 
interaction strength [5]. Moreover, the static properties of the 
system, such as the ground state, are also closely related to A. 
We can replace the operators fij with c numbers a% — y/nle 1 ^* 
in Eq. (1) to obtain a semiclassical Hamiltonian 
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(Nz 2 +N-2)- V z, (2) 



where <f> = <pi — <pR and z — (ul — nn)/N characterizing the 
imbalance. To minimize the energy, it is obvious that <f> should 
be zero. Then, for each Vq and A, we can find the position Zmi n 
of the local minimum of Eq. (2) by solving gj(Jf)|0=o = 0. 
For Vq < 10 -3 , at a critical A w 2, z m j n increases from 
abruptly, giving a hint of quantum transition of the ground 
state. 



III. FIDELITY SUSCEPTIBILITY 
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FIG. 1: (color online) The fidelity susceptibility x(A) for N = 800 
under tilt Vo = 10~ 10 , 10 -7 , 10~ 4 , 10 -3 and 10 _1 . One can see a 
clear transition from double peaks to single peak around Vo = 10 -5 . 



In this section, by diagonalizing the Hamiltonian (1), we 
use fidelity susceptibility to show the quantum transition of 
the ground state predicted by the semiclassical Hamiltonian 
(2). The ground state fidelity susceptibility is defined as 



y(A) = — lim 
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where \I/o(A) [vE f n (A)] is the ground (excited) state of %, 
E n (X) and [Eq(X)] is the corresponding ground (excited) en- 
ergy. Here we suppose that H has a non-degenerate ground 
state. In our system, ^ = — ^[hhip-L ~ 1) — !)]■ 

When Vq = 0, the ground states are degenerate for some val- 
ues of A. To break this degeneracy, we add an non-zero tilt 
and then we can use fidelity susceptibility to study the quan- 
tum transition of the ground state. 

In Fig. 1, we fix the particle number N to show the re- 
lation between xW an d Vo- We find that under a small Vq, 
x(A) has two peaks. With the increase of Vo, the position of 
the left peak does not change but the right peak moves left 
towards smaller A. When Vo is moderate (for example when 
Vo = 10~ 4 ), only one peak remains and continues moving 
left. When Vo is large enough (for example when Vo = 10 _1 ), 
no peak remains. 

The behavior of x(A) can be understood from the analy- 
sis of the ground state. The ground state can be expanded as 

l^o) = ^2k=o °k\k, N — k), where nj, = k and nji = N — k. 
Through studying the spectrum |cfe| 2 as a function of k, we 
can know the configuration of particles in the two wells. Fig. 
2 shows the spectra under different tilts Vo = 10~ 10 , 10~ 3 , 
and lfr 1 . Under V = 1(T 10 [Fig. 2(a)], when A < 2.06, 
|cfe| 2 is symmetric and has a peak at k = N/2, which is con- 
sistent with the knowledge that the ground state is a binomial 
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- k) at A = under 



V = 0. When 2.06 < A < 2.22, \c k \ z is still symmetric but 
has two peaks. This means the ground state becomes a cat-like 
state. When A > 2.22, \ck\ 2 only has one peak which moves 



towards k = N, meaning all particles tend to locate in one 
well and self-trapping occurs. Therefore, at small enough tilt, 
there are three phases, reflected by two peaks of fidelity sus- 
ceptibility. If we increase the tilt, the cat-like region becomes 
smaller and finally vanishes after the tilt is increased to an ap- 
propriate value, for example Vq = 10 -3 in Fig. 2(b). At this 
value of the tilt, the ground state will change from the bino- 
mial configuration to the self-trapping directly at some critical 
A, reflected by the single peak of x(A). For a very large tilt, 
\ck\ 2 is not symmetric in the whole region [Fig. 2(c)] and no 
ground state transition appears. 



So far our discussion is based on a fixed particle number 
N. Now we need to fix the tilt and enlarge N to see what 
happens in the large-iV limit. The height of peaks in both 
double-peak and single-peak region diverge with N exponen- 
tially as predicted by the fidelity susceptibility theory [Fig. 
3(c)]. In the double-peak region [Fig. 3(a)], when we increase 
N, the distance between the two peaks becomes smaller. A 
finite-size analysis of the positions A max of both peaks shows 



that |A max — 2| oc N~ d " with d p m 0.6799 for the left peak 
and d p pa 0.738 for the right peak. Therefore, A max of both 
peaks will tend to 2 when N — > oo, being consistent with 
the prediction of Eq. (2). Similarly, in the single-peak re- 
gion [Fig. 3(b)], |A max — 2| is also proportional to N~ dp with 
d p Ri 0.8941. We can know from these results that in the 
thermodynamic limit there is only one direct quantum phase 
transition from the binomial state to the self-trapping state. 
The double-transition from the binomial state to the cat state 
then to the self-trapping state under a small tilt and for mod- 
erate N is actually a crossover. Our results confirm that the 
fidelity susceptibility is useful for detecting not only quantum 
phase transitions but also crossovers [19]. 




FIG. 3: (color online) (a) The fidelity susceptibility x(A) for different particle numbers N =800, 1000 and 1200 under a small tilt Vo = 10 -10 . 
The inset shows the finite-size scaling analysis of the position A max of each peak. One can see that for both peaks, ln|A max — 2| oc IniV 
with a negative slope [green dots for the left peak (peak 1) and magenta dots for the right peak (peak2)]. Therefore in the large- N limit, both 
peaks are at the same position A max = 2. (b) The fidelity susceptibility xW f° r different particle numbers N =800, 1000 and 1200 under a 
moderate tilt Vo = 10 -3 . The inset shows the finite-size scaling analysis of the position A max of the single peak. Here one can also find that 
In | A max — 2| oc In N with a negative slope and in the large- N limit A max = 2. (c) The finite-size scaling analysis of x(A max ) under small 
and moderate tilt, which diverges exponentially with N. 



IV. QUANTUM DISCORD 

It is usually stated that the quantum transition of the ground 
state can be indicated by some quantum information quantity, 
such as entanglement of the ground state. Here we revisit this 
problem by studying the correlation in the transition region 
of our system. Before we discuss this, we first introduce two 
quantities we use to describe the correlation: the total correla- 
tion and quantum discord. 

Suppose that we have a system AB composed by two sub- 
systems A and B. Then we can use three density matrices 
Pab, Pa and p B to describe the states of the whole system and 
the two subsystems respectively, where p A m) = Tib(A)Pab- 
The joint entropy of the whole system is defined as the von 
Neumann entropy of p AB : S(p AB ) = -Tr(p AB \n p AB ). 
Similarly, we can also calculate the von Neumann entropy 
S{pa(b)) of Pa(b)- If Pab is a pure state, S(p A(B) ) = 
— Tt(pa(b) ^Pa(b)) is called as entanglement entropy and 
used to quantify the quantum entanglement between A and 
B. However, if p AB is a mixed state (like the two-body re- 



duced density matrix in our system), S(pa(b)) is not a good 
measure of entanglement. The total correlation (the quantum 
mutual information) between A and B is given by 1 AB — 
S{p A ) - S{p A \p B ), where S(p A \p B ) = S{p AB ) - S(p B ). 
Generally speaking, a bipartite quantum state p AB has both 
classical and quantum characteristics. So we can divide the 
the total correlation 1 AB into two parts: the quantum part and 
the classical part. 

The classical part is defined as the maximum information 
about one subsystem that can be obtained by performing mea- 
surements on the other subsystem. Let us consider a mea- 
surement performed only on subsystem B. This measurement 
can be described by a complete set of projectors {-M^} where 
Mk > and ^2 k Mu =Ib- The state of system AB after the 
application of M k becomes p k AB — ® Mk)pAB$A ® 

Mk) with p k = Tr[0U ® M k )pAB^A ® Mk)]. Accord- 
ing to the definition, the classical correlation can be ob- 
tained as C AB = S(p A ) - ™n{M k }J2kPk s (PA)> where 
Pa = T^bPab an d p\ = Tr B p AB . The quantum compo- 
nent of the correlation between two systems can be regarded 
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as the difference between the total correlation and the classi- 
cal correlation. This quantity is what we call quantum discord 
T>ab = 1-ab —Cab- It is interesting that some non-entangled 
states can also have nonzero quantum discord [14-18], reveal- 
ing that quantum discord is more suitable to capture the quan- 
tumness of correlation than entanglement. For pure states, the 
discord reduces exactly to the entanglement entropy. 

Now we consider quantum discord and the total correlation 
between two particles in our double-well system (See Ref. 
[20] for the discussion of the entanglement between identical 
particles). Because all particles are identical qubits, we have 
Tab — 2Si — S 2 , where Si(S 2 ) is the von Neumann entropy 
of one- (two-) particle reduced density matrix. Here the one- 
particle and two-particle reduced density matrices are defined 

as (Pi)ij = w( a ] a t) and (P2)ij,kl = Ni N-i) ( a Wi a j a ^) re- 
spectively, where i, j, k, I £ {L, R} and the average is made 
under the ground state. For qubits, each complete set of 
projectors contains two elements labeled by two parameters 
9 £ [0, 7r] and (p £ [0, 2tt], such that 



M 1 (e, v ) = 

M 2 (0,<p) = 
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One can note that actually we can write A4i(9,tp) 
\*i(e,tp))($i(0,tp)\ (i = 1,2) with 



cos I 
sin fe^ 



and \® 2 (8,<P)) = 
quantum discord can be expressed as 



sin : 



vp 



Then the 



T>AB = Si — S 2 
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min > 

k—1 



p k {9,ip)S{p k 2 {9,<p)], 



where Pk (0,<p) = Tr[(I <g> M k (0, <p))p 2 (I ® M h (0, <p))] and 
(*&<P) = ^- ) {l®M k {6,v)) P 2{I®M k {6,v)). Al- 
though in some cases one can obtain an analytical expression 
for T>ab [17], we have to do a numerical calculation here. We 
divide the domains of 9 ([0, ir]) and tf> ([0, 2tt]) into 100 equal 
intervals respectively and search the minimization. 

Similar with what we did in Sec. Ill, we first fix the particle 
number N to study the relation between correlations and Vb 
(Fig. 4). We find that T>ab and Tab have similar behaviors. 
Under a small tilt Vq = 10~ 10 , each correlation in the transi- 
tion region is remarkably larger than that out of the transition 
region, which means this transition cannot be described by the 
mean field theory for this system size. However, for each cor- 
relation we only observe one peak, whose position is near the 
position of the right peak of x(A) [Fig. 3(a)]. If we increase 
the tilt to Vb = 10~ 3 , both T>ab and Tab become smaller. 
However, a peak still exists, whose position is near the posi- 
tion of the single peak of x(A) [Fig. 3(b)]. If the tilt is further 
increased, both T>ab and Tab are almost zero, meaning the 
ground state is almost a product state for any A. 

Then, we want to know whether we can have a non-zero 
correlation in the large-TV limit. After fixing Vb, we find both 
T>ab and Tab decrease with N [Fig. 5(a) and (b)]. Through 
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FIG. 4: (color online) The quantum discord T>ab and total correla- 
tion Tab for N = 800 under different tilts Vo = 10" 10 , 10" 3 and 
10 _1 . One can see that with the increase of Vo, both T>ab and Tab 
decrease. 



a finite-size analysis for the tilt Vb = 10~ 10 , we find the peak 
values of the correlations decay exponentially as N~ dc (d c ps 
0.6661 for quantum discord and d c w 0.7426 for the total 
correlation). This means in the N —> 00 limit, there will not 
be correlation in the transition region and this transition can be 
described by the mean field theory. For Vb = 10~ 3 , a similar 
conclusion can also be obtained. 



V. SUMMARY 

In this paper, we analyze the quantum transition of the 
ground state for the single level Bose-Hubbard model in a 
double-well system with an arbitrary tilt. A semiclassical 
Hamiltonian predicts that for a not too large tilt, this transition 
happens at NU/J = 2. We use fidelity susceptibility x(A) to 
identify this transition. We find that for a small tilt, x(A) has 
two peaks which are at the same position in the N — > 00 limit. 
One peak corresponds to the transition from a binomial state 
to a cat-like state and the other peak corresponds to the transi- 
tion from a cat-like state to the self-trapping state. While for 
a moderate tilt, only one peak of x(A) is observed, which cor- 
responds to the direct transition from a binomial state to self- 
trapping. For a large tilt, no transition is observed in x(A). 

We also use two quantities describing correlation, quantum 
discord and the total correlation, to indicate the ground state 
transition. For a finite system size, each correlation has a peak 
in the transition region (either for a small tilt or a moderate 
tilt), meaning the transition cannot be described by the mean 
field theory. However, by doing a finite-size analysis, we find 
that in the N — > 00 limit, both correlations decay exponen- 
tially with particle number to zero. It is an interesting gen- 
eralization to check the behavior of quantum correlation with 
system size in other systems where a ground state transition 
that cannot be described by the mean-field theory exists, such 
as rotating BEC. 
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FIG. 5: (a) The quantum discord T>ab and the total correlation Tab for different particle number TV = 800, 1000 and 1200 under a small tilt 
Vb = 10" 10 . One can find that with the increase of TV, both Dab and Tab decrease, (b) The quantum discord T>ab and the total correlation 
Tab for different particle number TV = 800, 1000 and 1200 under a moderate tilt Vb = 10~ 3 . One can also find that with the increase of TV, 
both Dab and Tab decrease, (c) The finite-size scaling analysis of the peak values of both T>ab and Tab under the small tilt Vq = 10~ 10 
from TV = 3000 to TV = 9000. One can see that lnD^jf (InX^jf ) oc In TV with a negative slope, where (Tab*) is the peak value of 

Dab (Tab)- This means both of them decay to zero in large-TV limit. A similar analysis (not shown here) demonstrates that under a moderate 
Vb = 10 -3 , the peak values of both Vab andl^s also decay to zero in large-TV limit. 
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